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ABSTRACT 

We report findings from tlie first set of data in a current survey to estab- 
lisfi conclusively whether jets from young stars rotate. We observed the bi-polar 
jets from the T Tauri stars TH28 and RW Aur, and the blue-shifted jet from T 
Tauri star LkHa321, using the Hubble Space Telescope Imaging Spectrograph 
(HST/STIS). Forbidden emission lines (FELs) show distinct and systematic ve- 
locity asymmetries of 10 - 25 (± 5) km s^^ at a distance of 0".3 from the source, 
representing a (projected) distance of ~ 40 AU along the jet in the case of RW 
Aur, ^ 50 AU for TH28, and 165 AU in the case of LkHa321. These velocity 
asymmetries are interpreted as rotation in the initial portion of the jet where 
it is accelerated and collimated. For the bi-polar jets, both lobes appear to 
rotate in the same direction. Values obtained were in agreement with the pre- 
dictions of MHD disk-wind models (Bacciotti et al. 2002; Anderson et al. 2003; 
Dougados et al. 2003; Pesenti et al. 2003). Finally, we determine, from derived 
toroidal and poloidal velocities, values for the distance from the central axis of 
the footpoint for the jet's low velocity component of ~ 0.5 - 2 AU, consistent 
with the models of magneto-centrifugal launching (Anderson et al. 2003). 
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1. Introduction 

A key question in star formation research concerns the mechanisms behind the launch 
of jets from young stars. These jets are beheved to play an important role in the removal of 
excess angular momentum from the system, thus allowing accretion of matter onto the star 
up to its final mass. It is generally acknowledged that magneto-centrifugal forces are respon- 
sible for jet launching. In particular in the so-called 'disk-wind' model (e.g. Ferreira 1997; 
Konigl & Pudritz 2000) the jet is launched from the disk surface within a few AU from 
the star; while in the 'X-wind' model (Shu et al. 2000) the base of the flow is located at a 
few stellar radii from the source. To date, resolution constraints on observations have im- 
peded progress in validating the magneto-centrifugal mechanism, since jet launching occurs 
on small scales (i.e. less than 20 AU from the star); moreover, infall and outflow kinematics 
are complex and confused close to the source, that is often heavily embedded. Recently, 
however, interesting results have been obtained from observations of jets from more evolved, 
less embedded T Tauri stars (TTSs) for which the jet can be optically traced back to its 
origin. 

Observational backing for canonical models would require, for example, proof of ro- 
tation around the symmetry axis, close to the base where the jet is launched. The first 
hint of jet rotation was reported for the HH 212 system (Davis et al. 2000). However 
the knots were located at 2x10^ - 10^ AU from the jet source. Independently, asym- 
metries in velocities within the first 110 AU of the outflow from the T Tauri star DG 
Tau were found (Bacciotti et al. 2002), indicative of rotation. These results were obtained 
through an analysis of high angular resolution spectra taken with the Hubble Space Tele- 
scope Imaging Spectrograph (HST/STIS), aimed at probing the acceleration and collima- 
tion region of a stellar jet. Further confirmation of the rotation hypothesis came from 
Owens Valley Radio Observatory (OVRO) observations, which report the sense of rota- 
tion of the disk of DG Tau to be the same as that of the jet (Testi et al. 2002). More- 
over, the derived toroidal velocities in the observed portion of the jet were seen to be 
in agreement with the predictions of the magneto-centrifugal models, and indeed they 
can be used to find the location on the disk plane of the launching point of the wind 
(Bacciotti et al. 2002; Anderson et al. 2003; Dougados et al. 2003; Pesenti et al. 2003). 

These results motivated us to conduct an optical STIS survey to estabhsh conclusively 
whether jets from young stars rotate. We report findings for another three sources (of eight 
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in the survey sample) for which the data have already been acquired, i.e. the bi-polar jets 
from the TTSs TH28 and RW Aur, and the blue-shifted jet from LkHQ;321. 

2. Observations 

Observations were made of the jets associated with TH28, LkHQ;321 and RW Aur using 
HST/STIS on June 22, August 20, and October 3, 2002 respectively (proposal ID 9435). An 
aperture of 52 x 0.1 arcsec^ was used with the G750M grating which gave a spectral samphng 
of 0.554 A pixel"^, corresponding to a radial velocity of ~ 25 km for the wavelength range 
covered. The angular sampling was 0".05 pixcl^^. An acquisition of the stellar peak intensity 
prior to science observations allowed the slit to be centered accurately on the star before being 
offset to a position perpendicular to the jet axis at 0".3 from the source. We assume incli- 
nation angles (w.r.t. the plane of the sky) of 10° for TH28 (Krautter 1986), 44° for RW Aur 
(Lopez- Martin et al. 2003), and 45° for LkHQ;321 (arbitrary, since unknown, although it may 
have a large inchnation angle from spectro-astrometric measurements (Whelan et al. 2003)). 
The offset of 0".3 then represents a deprojected distance of ^ 51, 195 and 233 AU along the 
jet for TH28, RW Aur, and LkHa321, respectively. The only exception was in the case of the 
RW Aur blue-shifted jet lobe where the slit was placed at 0".2 due to lack of line emission 
at 0".3 (Woitas et al. 2002). Spectra were obtained of the blue- and red-shifted lobes, using 
exposure times of 2200 and 2700 seconds respectively, for the bi-polar jets associated with 
TH28 and RW Aur. In the case of LkHQ;321, only the blue-shifted lobe was detectable and 
so two spectra of this lobe, with exposure times of 2200 and 2700 seconds, were obtained and 
summed to increase the signal-to-noise of the faint emission lines. In total, this yielded five 
spectra, in the transverse direction at the base of the jets, which included Ha and the forbid- 
den emission lines (FELs) [OI]AA6300,6363, [NII]AA6548,6583, [SII]AA6716,6731. Data were 
calibrated through the standard HST pipeline, subtraction of the reflected stellar continuum 
was performed, and hot /dark pixels were removed. 
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Table 1: Sources studied in tiiis paper. Where known, the table also lists the mass of the star, tiie jet inclination angle 
ijet with respect to the plane of the sky, the mass accretion rate through the disk (Mdjajj), and the fluxes of mass (Mjet) and 
poloidal momentum (Pjet) in the jet. References: 1 - Bacciotti & Eisloffel 1999; 2 - Krautter 1986; 3 - Woitas et al. 2002; 4 - 
Lopez-Martm et al. 2003; 5 - Mundt & EisloflPel 1998. 
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3. Results 

In order to determine whether rotation is present in the jet channel, the difference in 
velocities on either side of the jet axis was analysed. Position- velocity contour plots for a 
sample of emission lines are shown in Figure 1, top panels. All radial velocities are systemic, 
i.e. they are quoted with respect to the mean heliocentric velocity of the star, which has 
been measured from photospheric lines to be +5 km s~^ for TH28 (Graham et al. 1988), 
+23 km s-^ for RW Aur (Woitas et al. 2003) and -7 km s"^ for LkHa321 (E. Whelan, priv. 
comm.). The lower order contours trace the outer jet channel where the jet is not so well 
coUimated and where the lower velocities lie. If rotation is present, there will be a difference 
in radial velocities between the two sides of the jet. This difference will be evident graphically 
as a skew in the contours of the transverse position-velocity diagram. Such a skew is indeed 
observed in all three cases in the outer contour lines. This suggests the presence of rotation 
in at least the low velocity component (LVC) of the flow near the outer borders of the 
jet channel. The high velocity component (HVC), which is located much closer to the jet 
axis and gives rise to the emission peak (Bacciotti et al. 2000), appears not to be spatially 
resolved in our spectra. For this reason we cannot resolve any velocity difference in the 
HVC between the the two sides of the flow, i.e. we cannot detect rotation for this velocity 
component. 

To give a quantitative estimate of the observed velocity shifts, it was flrstly necessary 
to ensure that we measure velocity offsets at equal distances on either side of the jet axis. 
To this purpose we assumed that the peak of the (HVC) emission traces the position of 
the axis, and measured its distance from the nominal center of the slit, with a Gaussian 
fit along the cross-dispersion direction. We then shifted the line emission re-centering the 
HVC peak on the nominal zero arcsecond position. In all cases a small offset (< 0.4 pixels) 
was required. Table 2. The offsets have different sign and magnitude for the three targets, 
indicating that this is not an instrumental systematic effect. Instead, pixel shift values 
for emission lines in opposite jet lobes are consistent with a maximum misalignment < 2 
degrees with respect to the perpendicular of the actual jet axis at sub-arcsecond scales. This 
resulted in a displacement of the peak intensity of up to 0".02, a quantity which would 
produce a marginal asymmetry in radial velocity estimated to be at most 20% if uncorrected 
(Bacciotti et al. 2002). This error was, however, avoided as the emission was brought on-axis 
prior to analysis. 

The peak intensities of pixel rows on either side of the central row were then compared 
for velocity differences. Figure 1, bottom panels. The intensity profiles of each pair of pixel 
rows symmetric about the jet axis was plotted (e.g. ±0".l, or as indicated in each box). The 
single curve at the bottom in each case is the intensity of the central on-axis pixel row. Two 
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Fig. 1. — Top row: Selected position- velocity contour plots of emission lines from the studied 
jets. The skew in lower order contours is indicative of rotation in the outer jet channel, while 
the high velocity component remains unresolved. Contour values in units of erg cm~^sec~^ 
A"^ arcsec"^ in the panels are - TH28: from 7.0 x 10"^^ to 6.7x10"^^ with intervals 3.0 x 10"^^; 
RW Aur: from 7.2x10-^^ to 2.3x10-^^ log scale interval of 2^/^; LkHa321: from 1.0x10"^^ 
to 8.2x10"^^ with intervals of 4.0x10"^^. Panels are corrected for the heliocentric velocity 
of the star, that is for +5, +23 and -7 kms"-*^ in the three cases, respectively. Bottom 
row: Normalized intensity profiles along horizontal cuts in the above panels, symmetrically 
opposed with respect to the jet axis. 
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methods were used in velocity measurements: a cross-correlation technique, which analyses 
the overall displacement of lines and is independent of the shape of the line profile; and a 
Gaussian fitting technique, which acts as a suitable check given the simple shape of the line 
profile in most cases. Specifically, each pair of pixel rows, mirrored in distance from the jet 
axis, was cross-correlated, and Gaussian fits for each pair of rows were also compared. The 
outcomes of the two methods are consistent, showing clear radial velocity differences of 10 
to 25 (± 5) km s"-*^ for opposing jet edges. Results are listed in Table 3, in which we report 
the measured radial velocity differences in the direction of the oriented slit, corresponding 
to a direction specified in the first column of the Table. In a small number of cases, fiagged 
in the table with a symbol, the emission had to be filtered out of a background which was 
causing velocity measurements to be artifically changed due to either: relatively strong HVC 
presence; oversubtraction of the background; or low signal-to-noise at crucial positions. The 
overall results are illustrated graphically in Figures 2 to 4. 

Finally, the radial velocity profile across each jet is shown in Figure 5. In the red-shifted 
lobes of TH28 and RW Aur, the on-axis radial velocity is clearly the highest and the value 
reduces as the edges of the jet are approached. The blue-shifted lobes are not so clear as the 
emission is fainter and often more scattered. Data reduction in these cases often required 
isolation of the lower velocity component from the HVC (as marked in Table 3) to identify 
elements of rotation. However radial velocities in this figure include the HVC, and so rotation 
is not apparent in all cases. 

From the results of this spectral analysis, combined with the inclination angles, we find 
poloidal velocities for the RW Aur jet of 144 - 227 km s^^ in the red lobe and 245 - 288 
km s~^ in the blue lobe - this velocity asymmetry is well known from previous observations, 
(Woitas et al. 2002); for TH28 jet of 115 - 288 km s"^ in the red lobe and 230 - 374 km s'^ in 
the blue lobe - an asymmetry which was also previously recorded, (Graham et al. 1988); and 
for LkHQ;321 of 540 - 550 km s~^. Toroidal velocities derived from the outer positional radial 
velocity shifts, being less affected by projection effects (see Section 4), are in the ranges of: 
7 to 17 km s"^ for both lobes of RW Aur; 5 to 13 km s"^ for the red lobe of TH28; 4 to 8 
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Table 2: Pixel shifts applied to emission lines in each jet lobe along the cross-dispersion direction, to re-center the HVC 
emission peak, assumed coincident with the jst axis, on the nominal zero arcsecond position. 
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Fig. 2. — Left: schematic drawing of the observing mode for the bhie-shifted jet from 
LkHQ!321. The derived sense of rotation is illustrated by the circular arrow around the 
symmetry axis. Top right: Symbols adopted for the various lines. Bottom right: Summary 
of radial velocity asymmetries measured at the base of the blue-shifted lobe, at 0".3 from 
the source, in the direction of the oriented sht, i.e. North West - South East (also specified 
by the label of the y-axis). 
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Fig. 3. — Same as Fig. 2, for the bi-polar jet from TH28. Both shts were located at 0."3 
from the source. 
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Fig. 4. — Same as Fig. 2, for the bi-polar jet from RW Aur. In the blue-shifted jet lobe the 
spectrum was taken at 0".2 from the star. 
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Tabic 3: Radial velocity differences, AVrad- across the jet 0".3 from the source (0".2 for the RW Aur blue lobe), measured 
using both the eross-eorrelation teelmique and single Gaussian fitting, along the direction specified in the first column. The 
symbol * marks data points which have been filtered out of a low S/N environment, or from a strong spurious contribution 
from the HVC emission wings. Where dots appear in the table, the emission was either shifted off the CCD, or was too faint, 
small or scattered to decipher. The accuracy reached with the data analysis is approximately ±5 km s~^. 
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km s-^ for the blue lobe of TH28; and 4 to 9 km s"^ for the blue lobe of LkHQ;321. (Note 
that, in the ranges given above, the higher poloidal velocities lie at distances closest to the 
rotation axis, and these correspond to the lower toroidal and radial velocities.) 

4. Discussion 

Assuming that emission from the jet is axially symmetric, we interpret our findings 
of velocity differences between the two sides of the jet flow. Table 3 and Figures 2 to 4, 

as evidence for rotation at the base of the jet. Most importantly, for TH28 and RW Aur, 
the red and blue jet lobes were found to rotate in the same direction. This implies that 
the helicity in the red and blue lobes (i.e. the handedness of toroidal with respect to the 
poloidal velocity) is opposite in opposite directions, a result predicted by MHD models where 
the ambient field is wrapped around due to disk rotation. 

Before looking at each target individually, there are a few general comments to be 
made about the results in Table 3. Firstly, the velocity difference measurements close to 
the jet axis, i.e. at the 0".05 position, are noticeably smaller than points further away, and 
where the emission is detectable at 0".2 and beyond, velocity differences are higher than at 
intermediate distances. This effect may at first appear in contradiction with the notion that 
the central portions of the jet should rotate faster. A detailed comparison with disk-wind 
model predictions (Pesenti et al. 2003, and Dougados et al. 2003) shows, however, that the 
apparent decrease of the observed velocity difference toward the jet axis is likely to be due to 
projection and beam smearing effects. Since the emission is optically thin we see, along the 
line of sight, the sum of the contributions emitted from regions in the jet that rotate with 
different toroidal velocities. This causes a reduction of the observed velocity shift. Such an 
effect is more important for regions closer to the axis, while the values measured at the outer 
jet borders are less contaminated and so are in better agreement with theoretically predicted 
toroidal velocities. In other words this effect does not reflect a true kinematic feature, but 
is expected on the basis of MHD acceleration models when combined with our observational 
mode. Secondly, it should be noted that the size of velocity differences in different emission 
lines does not represent scattering around an average value but rather is due to the fact that 
emission has its origin at different positions along the line of sight. And lastly, the [Nil] lines 
show higher velocity differences than other emission lines, a result that illustrates how they 
trace the central more collimated higher velocity region of the flow (see Bacciotti et al. 2000, 
Lopez-Martin et al. 2003, Pesenti et al. 2003). 

The results for the TH28 red-shifted jet lobe are clearest. Velocity differences alomg 
the SW - NE direction are positive, with a few exceptions mainly in the [SIX] lines. Also, 
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Fig. 5. — Radial velocity profile across each jet in various emission lines. The TH28 blue- 
shifted [01] A6300 line is not included as the profile of this emission line is very wide and 
did not allow measurement of velocity peaks with Gaussian fitting. In this case the rota- 
tional velocities were derived from cross-correlation routines alone, resulting only in velocity 
difference measurements. 
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the data close to the jet axis are not well resolved as explained above. The stronger [01] 
and [Nil] lines have values of 5 and 6 km s~^ at 0".05 from the axis compared to 10 to 20 
km s~^ further from the axis, while the outer jet channel seems to have higher radial velocity 
differences of about 24 km s^^. Values in the blue-shifted jet lobe are less clear, but [Nil] and 
[01] emission give positive differences consistent with the red-shifted lobe. All other usable 
data points fall within the error bars about zero. Globally, we find evidence that both lobes 
of the jet rotate in a clockwise direction, looking down the blue-shifted jet lobe towards the 
source, with a measured radial velocity difference of 10 to 25 km s~^. 

The RW Aur jet also shows clear indications of rotational velocities. Exceptions include 
the [S11]A6731 values in the red-shifted lobe. Also points close to the jet axis show smaller 
radial velocity differences due to strong unresolved HVC emission. However, the [01] hues 
give clear results with higher radial velocity difference evident at 0".2 from the jet axis. 
The blue-shifted lobe is less definite, but nevertheless velocity differences outside the error 
bars are positive in line with the red-shifted lobe. Overall, results show an anti-clockwise 
rotation looking down the blue-shifted jet lobe towards the source, again with radial velocity 
differences of 10 to 25 km s^^. These findings are consistent in magnitude and direction 
with results of similar research on the RW Aur jet (Woitas et al. 2003) in which rotational 
velocities in the same direction of 10 to 20 km s~^ have been observed, in the form of velocity 
differences between the borders of the flow. For this study HST/STIS was also used, but 
the spectra were taken in a set of positions across the jet with the slit direction parallel to 
the jet axis. 

In the case of LkHa321, which is located at 550 pc (more than three times the distance 
of the other targets) the emission lines were very faint despite having combined two spectra 
to increase the signal-to-noise ratio. Nevertheless, velocity differences in the range of -5 to 
-30 km s~^ have been measured, although evidence of rotation is weaker in this case. The 
sense of rotation for LkHQ;321 is measured as anti-clockwise looking down the blue-shifted 
lobe towards the star. In this jet we are able to measure rotation at higher velocities because 
of a number of factors: LkHQ!321 is farther away meaning we are looking further along the 
jet (233 AU from the source, projected distance) to a point where it has widened and so we 
can resolve higher velocities; a large inclination angle (which appears to be the case from 
spectro-astrometric measurements (Whelan et al. 2003)) would mean that the poloidal and 
radial velocities approach each other; and finally this is a larger mass T Tauri star, since 
it has spectral type Gl (Chavarria-K. et al. 1981), implying higher velocities in the outer 
resolvable regions of the jet. 

It is reassuring that we obtain negative shifts in the case of LkHQ;321, as opposed to 
the positive shifts observed for the other two sources because, as we are comparing the same 
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rows on the CCD detector in all cases, it means we are not measuring an instrumental effect. 
Conceiveably, a slight misalignment of the slit with respect to the transverse direction of the 
jet (i.e inaccurate position angle used in the pointing of the instrument) may produce an 
effect similar to rotation. In this case, the position of the real jet axis will then be shifted 
with respect to the nominal zero arcsec position. We have shifted the spectral image back to 
the zero arcsec position as previously discussed (Section 3), but the angle subtended remains 
a problem. It could produce a rotation signature even in a non-rotating jet, since the HVC 
and LVC are at different spatial locations on the CCD with respect to the zero arcsec row 
of pixels, implying that we are not probing symmetric regions of the jet with respect to its 
real axis. Given such a misalignment has occurred, the extent of the contamination does 
not have, however, a dramatic affect on our results, as can be seen from the fact that the 
pixel shift requirements for TH28 and LKHa321 (Table 2) are in the same direction but 
the sense of rotation of their jets is opposite. For RW Aur, the same PA was used here 
for slit positioning as in a previous study (Woitas et al. 2003), where it was found that the 
magnitude of the false rotation signature contamination due to incorrect PA was, at most, 1 
- 5 km s""^. However, the sense of the false signature is in fact opposite in direction to that of 
the jet's rotation, and so the values we have measured are actually lower limits. Apart from 
this, the only other obvious effect which could produce a contour skew mimicking rotation is 
asymmetrical interaction with the local environment on either side of the propagating jet, e.g. 
asymmetrical mass entrainment leading to asymmetrical poloidal velocities. However such 
mimicking is unlikely because (apart from the fact that asymmetrical entrainment should 
also produce enhanced emission at one border of the flow, which is not seen in our spectra) 
we see the same asymmetry in both the red and blue shifted jet lobes where present, and the 
magnitudes of the deduced toroidal velocity differences are in the range predicted by theory 
(as discussed below). 

Overall, our observations are in line with the observations of the jet from the T Tauri 
star DG Tau (Bacciotti et al. 2002). In that case, by using simple and general relationships 
governing the physics of magnetically launched disk winds together with an observationally 
based estimate of the ratio, R, between the mass flux in the jet and the mass flux accreted 
through the disk {R ~ 0.1), it was demonstrated that the observed velocity differences were 
in the expected range, (Bacciotti et al. 2002: Anderson et al. 2003; Dougados et al. 2003). 
These values compare well with our results, which therefore support the magneto-centrifugal 
scenario. Furthermore, our toroidal and poloidal velocities have the same ratio as theoretical 
predictions (Vlahakis et al. 2000), and we can use these velocities, v^^oo and Vp^oo measured 
at a distance, ru^, from the rotation axis, to obtain values for the wind-launch region 
in terms of distance from the rotation axis along the disk plane, Wq. Using Equation 1 
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(Anderson et al. 2003), 



2/3 / Vrhr^ ^2/3 

1/3 



X 



and assuming all three sources are of mass ~ IMq, which is a reasonable approximation 
given the weak dependence on M^, we obtain values for wq as shown in Table 4. We have 
chosen measurements at 0".l and 0".2 from the jet axis as hmits of a suitable range, being 
less contaminated by projection effects. Where emission was faint, values corresponding 
to 0".05 were used, but it should be noted that these are less accurate. Also, the values 
at 0".l are less precise than those at 0".2, when both are present, since increased projec- 
tion effects close to the jet axis tend to reduce the line of sight averaged v^^^o, and hence 
the resulting value of Wq (Pcscnti et al 2003). Considering the red-shifted jet-lobes which 
both have very strong signal-to-noise in [01] and [Nil] emission lines, we calculate a wind- 
launch region spanning 0.3 to 1.6 AU, in keeping with previous estimates for DG Tau of ^ 
1.8 AU (Bacciotti et al. 2002) and ~ 0.3 to 4 AU (Anderson et al. 2003). The observational 
evidence presented here supports the idea that disk winds are launched, via the magneto- 
centrifugal mechanism (e.g. Konigl & Pudritz 2000), at radii within a few AU of the star. 
Determination of the most appropriate model however will have to await higher spatial and 
spectral resolution observations in the future. 



5. Conclusions 

The jets from the three young stars observed, TH28, RW Aur and LkHQ;321, show 
distinct and systematic radial velocity asymmetries in opposing positions with respect to 



Star 


Jet lobe 








'^radial 








arcsec 


AU 


km s"'^ 


km 


AU 


TH28 


Red-shifted 


0.1 - 0.2 


17 - 34 


10 - 25 


30 - 20 


0.3 - 1.6 




Blue-shifted 


0.1 - 0.2 


17 - 34 


8 - 15 


50 - 40 


0.1 - 0.4 


RW Aur 


Red-shifted 


0.1 - 0.2 


14 - 28 


10 - 25 


105 - 100 


0.4 - 1.3 




Blue-shifted 


0.05 - 0.1 


7 - 14 


10 - 25 


180 - 170 


0.1* - 0.4 


LkHa321 


Blue-shifted 


0.05 - 0.1 


27.5 - 55 


5 - 12 


390 - 380 


0.1* - 0.2 



Table 4: The range for the launch point of the disk wind, Wq, for our five targets, calculated 
using the method described in Anderson et al (2003). *Where emission is faint, values for 
the 0".05 position were used, but these should be considered less accurate due to resolution 
constraints. 
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the jet axis, at 0".2 - 0".3 from the source. Although the on-axis higher velocity component 
of the jet remains unresolved, radial velocity differences in the lower velocity component 
located in the outer jet channel are found to be on the order of 10 to 25 (±5) km s~^. 
For the bi-polar jets from TH28 and RW Aur, the velocity differences have the same sign 
in both lobes. We interpret these radial velocity asymmetries as rotation signatures in the 
region where the jet has been collimated but has not yet manifestly interacted with the 
environment. Therefore the sense of rotation of the jets, looking down the blue-shifted lobe 
towards the star, is clockwise for TH28, and anti-clockwise for RW Aur and LkHQ;321. 

Our findings arc reinforced in a number of ways: the velocity differences arc of the same 
magnitude as those measured in the similar DG Tau jet of 5 to 10 km s^^ (Bacciotti et al. 2002), 
which was shown to be in agreement with the predictions of MHD disk-wind models (Bacciotti et al. 2002; 
Anderson et al. 2003; Dougados et al. 2003; Pesenti et al. 2003); they are in hne with sim- 
ilar research on the RW Aur jet (Woitas et al. 2003) which yields rotational velocities of 10 
to 20 km s~^, with the same sense of rotation; and finally, they lead to values for the distance 
of the LVC footpoint from the central axis of 0.5 - 2 AU, consistent with the models of 
magneto-centrifugal launching (Anderson et al. 2003). 
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